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NOMENCLATURE 

a  characteristic  length,  half  width  of  contact  zone  or  length  of  slider 
Cj  specific  heat  of  region  i 

D  Deborah  number,  n0ur/(Ga) 

d  surface  displacement 

d  dimensionless  surface  displacement.  2rd/a2 

E  elastic  modulus 

e  effective  strain  rate,  see  Equation  (20) 

F  derivative  function,  dx'/d5 

f  frictional  force  per  unit  transverse  length 

f  dimensionless  frictional  force,  f/(ak) 

G  shear  modulus 

h  coating  thickness 

h*  constant  related  to  flow,  2Q/(u0h) 

h*  dimensionless  constant,  h*/h° 

h0  characteristic  coating  thickness,  minimum  coating  thickness  under  linear  slider 
h1  maximum  coating  thickness  under  linear  slider 
fi,  dimensionless  coating  thickness,  f^/hg 

Kj  thermal  conductivity  of  region  i 

k  shear  yield  stress 

k  dimensionless  stress,  k/k0 

k0  characteristic  value  of  shear  yield  stress 

p  pressure 

p  dimensionless  pressure,  p/p0 

p  dimensionless  pressure,  ph0/(ak) 

ph  maximum  Hertz  pressure 

p0  reference  pressure,  taken  as  ph  for  Hertz  contact 
Q  volumetric  flow  per  unit  transverse  length 

heat  flux  crossing  surface  of  region  1,2 
qt  dimensionless  heat  flux,  q1/(4rurSTc0) 

P|  Peclet  number  for  region  i,  Uja/Kj 

Rj  volumetric  rate  of  heat  generation  for  region  I 
r  cylinder  radius 

slide  to  roll  ratio,  Au/ur 

dimensionless  shear  rate  parameter,  jigUfS/fhtgg) 
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residence  time  in  contact  zone 
temperature 

dimensionless  temperature  rise,  (T  -  T0)/Tr 
surface  temperature  at  start  of  contact 
reference  temperature  rise,  ursai(f 
velocity  in  x  direction 
dimensionless  velocity,  u/u0 
rolling  speed 

characteristic  velocity  in  x  direction 
speed  of  region  1 ,2 
velocity  in  y  direction 
dimensionless  velocity,  va/(hQu0) 
load  per  unit  transverse  length 
dimensionless  load,  Wh0/(a2k) 
coordinate  in  along  surface 
coordinate  perpendicular  to  surface 
viscosity  -  pressure  coefficient 
dimensionless  coefficient,  aph 
limiting  shear  stress  -  pressure  coefficient 
dimensionless  coefficient,  Pph 
viscosity  -  temperature  coefficient 
dimensionless  coefficient,  yTr 
inclination  of  slider,  h,  -  h0 

temperature  rise  on  surface  1  at  center  of  contact  for  uniform  flux 

temperature  rise  in  coating 

slip  rate,  u1  -  u2 

slip  velocity  at  y  =  h 

slip  velocity  at  y  =  0 

dimensionless  parameter,  (h0/a)(k0/p0) 

dimensionless  parameter,  maximum  of  ft  and  (a/h0)2 

shear  rate 

dimensionless  coordinate,  y/a 

thermal  diffusivity  for  region  i 

limiting  shear  stress  -  temperature  coefficient 

dimensionless  coefficient,  A.Tr 
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viscosity  at  ambient  pressure  and  temperature  T0 
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|i  dimensionless  viscosity,  p./p.Q 

A  dimensionless  viscosity-speed  parameter,  pu0/(h0k) 

A*  critical  value  of  A  beyond  which  W  =  0 
v  Poisson's  ratio 

i  dimensionless  coordinate,  x/a 

Pi  density  of  region  i 

o  effective  stress,  see  Equation  (19) 

ox,oy  coordinate  stresses  in  x,y  directions 
ox,oy  dimensionless  coordinate  stresses,  ox/p0,  oy/p0 
az  dimensionless  stress  in  transverse  direction 
x  shear  stress 

i  dimensionless  shear  stress,  x/k 

x  dimensionless  shear  stress,  x/xc0 

x  average  shear  stress  in  contact  zone 
x  dimensionless  average  shear  stress,  x/Tco 
t'  shear  stress  without  regard  to  limiting  shear  stress 
x'  dimensionless  shear  stress,  x7xc0 

x'  time  derivative  of  t' 

xc  characteristic  shear  stress  function  used  in  rheological  model 
xc  dimensionless  shear  stress,  xc/xc0 

tcQ  characteristic  shear  stress  at  ambient  pressure  and  temperature  T0 

xm  limiting  shear  stress 

tm  dimensionless  limiting  shear  stress,  Tm/xc0 

xh  shear  stress  at  y  =  h 

xh  dimensionless  shear  stress,  -xh/k 

x0  shear  stress  at  y  =  0 

x0  dimensionless  shear  stress,  -x0/k 

xmax  maximum  of  xQ  and  xh 

<t>  function  characterizing  rate  dependence  in  rheological  model 
X  limiting  shear  stress  ratio,  Tm/xc,  presently  taken  near  1 
i|r  flux  partition  fraction,  see  Equation  (11) 

Subscripts 

1  upper  region  (see  Figure  1) 

2  lower  region  (see  Figure  1) 

3  coating 

i  region  variable,  i  =  1 ,2,3 

j  particular  point  in  discretized  Z,  grid,  used  in  finite  difference  solution 
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I.  INTRODUCTION 

Solid  lubricants  have  been  in  use  for  many  years  in  situations  where  extreme  temperature,  pressure  or 
environmental  conditions  prohibit  the  use  of  liquid  lubricants.  Present  requirements  for  high  temperature 
and  hostile  environment  applications  create  an  increasing  need  for  improved  solid  lubrication.  Although 
considerable  strides  have  been  made  in  the  development  and  implementation  of  solid  lubricants,  they  have 
been  somewhat  hampered  by  a  lack  of  quantitative  methods  for  assessing  the  performance  of  solid  films. 

This  study  is  intended  to  represent  a  step  towards  the  objective  of  developing  a  tractable  means  for 
predicting  the  performance  of  solid  lubricating  films  in  machinery  based  upon  property  data  obtained  under 
laboratory  conditions.  Performance  characteristics  under  consideration  include  surface  and  sub-surface 
temperature,  friction,  operating  coating  thickness  and  wear  rates.  These  quantities  may  in  turn  be  used  to 
establish  the  useful  life  and  applicability  of  a  solid  lubricant  for  a  given  application. 

Considerable  emphasis  is  placed  here  on  the  need  for  a  tractable  approach.  Although  an  ultimate  goal  of 
predicting  performance  of  solid  lubricants  based  upon  molecular  structure  should  be  pursued,  a  shorter  term 
approach  is  also  necessary  to  satisfy  current  needs  even  though  it  might  require  additional  property  data. 
This  approach  involves  treating  the  solid  lubricant  as  a  continuum  and  modeling  the  rheology  of  the  contact 
based  on  properties  of  both  the  lubricant  and  its  immediate  surroundings. 

This  investigation  deals  with  an  idealized  isotropic  lubricant  under  conditions  of  plane  strain.  Normal 
stresses  due  to  loading  are  assumed  in  general  to  be  high  compared  with  the  shear  yield  stress  of  the 
lubricant.  The  film  is  assumed  to  be  thick  compared  with  asperity  heights  but  thin  compared  with  contact 
lengths.  In  the  first  part  of  the  study,  as  described  in  the  next  section,  the  existence  of  such  a  film  is 
assumed  and  efficient  algorithms  are  generated  for  predicting  tractions  based  on  contact  geometry, 
operating  conditions  and  the  rheological  model  for  the  lubricant.  Rheological  models  include  a  linear  elastic 
component,  a  non-linear  rate  dependent  component  and  a  limiting  shear  stress.  Properties  are  allowed  tc 
be  dependent  on  both  temperature  and  pressure.  A  thermal  analysis  is  included  and  coupling  between  the 
effects  of  heat  generation  and  temperature  dependent  properties  is  treated.  Numerical  examples  are 
presented  which  show  the  influence  of  the  various  parameters  on  shear  stress  behavior.  The  second  part 
of  the  study,  described  in  Section  III,  deals  with  the  ability  to  form  and  sustain  the  films  assumed  above  and 
limitations  that  occur  when  the  slip  rates  are  high  and  is  intended  to  provide  guidance  regarding  the  range 
of  validity  of  both  the  present  and  previous  traction  models. 
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II.  INITIAL  TRACTION  ANALYSIS 

A  combined  thermal-elastic-plastic-creep  model  is  presented  here  for  estimating  tractive  forces  resulting  from 
shearing  of  solid  films.  Models  of  this  type  have  been  used  successfully  in  the  past  for  predicting  tractive 
forces  arising  in  concentrated  contacts  with  liquid  lubricants.  Under  the  extreme  pressures  characteristic 
of  these  contacts,  liquids  behave  similarly  to  solids  in  that  elastic  properties  become  significant  at  low  shear 
rates  and  plastic  yielding  occurs  at  higher  shear  rates. 

It  will  be  assumed  that  the  film  is  thick  compared  with  the  surface  asperities  on  the  mating  surface  and  the 
load  is  sufficiently  high  that  the  asperities  are  fully  imbedded  in  the  film.  A  nearly  hydrostatic  state  of  stress 
will  be  assumed  to  exist  in  the  contact  zone.  The  normal  loading  will  affect  the  size  of  the  contact  zone  for 
elastic  substrates  and  the  pressure  distribution  within  the  contact  zone.  It  is  assumed  however,  that  the  only 
dependence  of  shear  forces  on  the  normal  loadino  arises  from  pressure  dependence  of  mechanical 
properties. 

Based  on  the  above  assumptions  the  initial  traction  model  will  be  an  adaptation  of  that  used  for 
elastohydrodynamic  applications1,2.  The  extensions  put  forth  here  will  include  coupling  with  a  fairly  general 
thermal  analysis  and  the  ability  to  treat  pressure  distributions  other  that  Hertzian. 

The  development  of  the  traction  model  contained  herein,  will  proceed  with  a  formulation  of  the  mechanical 
problem,  a  formulation  of  the  thermal  problem,  a  description  of  the  method  of  solution  to  the  combined 
problem  followed  by  some  parametric  studies. 

A.  Formulation  of  Mechanical  Problem 

The  analysis  developed  thus  far  is  based  on  a  two  dimensional  geometric  model  depicted  in  Figure  1 .  A 
rolling  -  sliding  contact  is  depicted  with  surfaces  moving  at  speeds  u,  and  u2.  A  contact  zone  is  formed 
between  x  =  0  and  x  =  2a  where  the  coating  thickness  is  assumed  to  have  a  constant  value  h.  The  shear 
stress  x  acting  on  the  upper  surface,  will  be  in  the  direction  shown  when  u,  >  u2.  A  class  of  rheological 
models  in  the  form 
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will  be  examined.  In  the  above  equation  G  represents  a  shear  modulus,  xm  is  a  maximum  limiting  shear 
stress,  and  xc  and  p  are  constants  and  4>  is  a  function  describing  the  creep  characteristics  of  the  coating. 
If  <(>  is  chosen  so  that  its  limiting  form  at  small  values  of  its  argument  approaches  its  argument  then  n  would 
reduce  to  a  viscosity  at  low  shear  stress.  In  general  all  of  the  properties  in  the  above  model  can  be 
dependent  on  temperature  and  pressure,  however  elastic  properties  will  be  somewhat  less  sensitive  than 
creep  and  plastic  properties  and  will  be  taken  to  be  constant.  This  includes  both  the  shear  modulus  of  the 
coating,  G,  as  well  as  the  elastic  properties  of  the  substrate. 


The  variation  of  pressure  across  a  thin  film  will  be  negligible.  For  thin  metallic  films  the  thermal  conductivity 
will  be  sufficiently  high  so  that  the  variation  in  temperature  across  the  film  will  be  small  when  compared  with 
the  variation  along  the  film.  If  properties  are  taken  to  be  constant  across  the  film  (evaluated  at  either  the 
surface  temperature,  the  centerline  temperature  or  a  mean  temperature)  then  the  velocity  distribution  will 
vary  linearly  across  the  film  and  the  strain  rate  will  be  given  by  e  =  Au/h,  where  Au  =  u,  -  u2. 

If  the  rolling  speed  ur  =  (u,  +  u2)/2  is  large  compared  with  the  slip  rate  Au,  one  may  approximate  the 
shear  stress  derivative  with  x'  =  ur  dx'/dx  and  the  top  part  of  Equation  (1)  may  be  written  as 
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The  pressure  and  temperature  dependence  of  properties  xc,  xm,  and  p  may  be  expressed  as 


TC-  tc0Xc(P.T),  Tm  -  TmO*m(P<T)  .  p-p0p(p,f), 

where  xc0,  xm0,  and  n0  denote  properties  evaluated  at  the  upstream  surface  temperature,  T0,  and  ambient 
pressure.  The  quantities  p  and  f  are  dimensionless  forms  of  the  pressure  p,  and  temperature  T,  scaled  on 
suitable  reference  quantities  which  will  be  selected  later. 

One  may  now  introduce  the  dimensionless  variables  and  parameters 
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to  put  Equation  (1)  in  the  dimensionless  form 


-  S  -  4r  *(*'/?„) 


t  -  . 


(4) 


If  we  neglect  shearing  effects  in  the  inlet  (5  <  0)  then  the  boundary  condition  on  Equation  (4)  is 

t  -  0  @  5  -  0  .  (5) 

If  properties  were  independent  of  pressure  and  temperature,  or  if  the  pressure  and  temperature  were  known 
as  functions  of  position  then  Equations  (4)  and  (5)  could  be  solved  for  any  prescribed  function  <J>.  Two 
particular  functions  to  be  considered  initially  are  4>(C)  =  tanh_1(C)  (to  be  referred  to  as  the  arctan  model), 
which  was  used  by  Bair  and  Winer  to  characterize  traction  fluids  and  4>(0  =  C  .  the  Maxwell  model  which 
has  been  used  extensively  in  the  analysis  of  viscoelastic  behavior  of  materials. 

The  following  representations  will  be  implemented  to  characterize  effects  of  property  variations  with 
temperature  and  pressure: 


H  ~  Hoe 


op  -  y(T  -  Tq) 


-  X, 


cO 


1  +  Pp 


1  +  X(T  -  T0) 


•m 


(6) 


In  order  to  couple  Equation  (6)  with  Equation  (4)  the  relationships  must  be  provided  for  determining  the 
pressure  p,  and  the  temperature  T,  as  functions  of  position  x.  Two  types  of  pressure  distributions  that 
should  be  considered  are  a  uniform  distribution  (whose  relevance  will  be  discussed  in  Section  III)  and  an 
elliptical  distribution  corresponding  to  a  Hertzian  contact  based  on  elastic  properties  of  the  substrate.  This 
latter  distribution  which  is  illustrated  in  Figure  1  will  be  used  in  the  initial  analysis.  Thus 


P  -  P„  Vi  -  K  -  i)2 


(7) 


where  ph  is  the  maximum  Hertz  pressure  easily  determined  from  the  Hertz  formulas  for  line  contact  along 
with  the  half  width  of  the  contact,  a.  The  use  of  this  relationship  will  be  valid  when  a/h  >  >  1  and  the 
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stiffness  of  the  film  in  the  contact  zone  is  high  compared  those  for  the  substrates. 

B.  Formulation  of  Thermal  Problem 


The  amount  to  heat  dissipated  in  the  contact  will  depend  on  the  shear  stress  distribution,  hence  the 
temperature  distribution  cannot  be  prescribed  apriori.  A  form  of  the  energy  equation  must  therefore  be 
coupled  with  the  rheological  equations  given  above.  The  energy  equation  will  first  be  simplified  based  on 
high  Peclet  number  approximations. 

Since  body  1  moves  at  a  speed  uv  the  residence  time,  tr,  of  a  point  on  its  surface  in  the  contact  zone 
(0  <  x  <  2a)  would  be  tr  =  2a/ur  If  the  thermal  diffusivity  of  body  1  is  then  the  thermal  diffusion  depth 
into  body  1  would  be  Vk  1 1  r  and  the  ratio  of  the  length  of  the  contact  zone,  2a,  to  the  diffusion  depth 
would  be  ^2Pt ,  where  P,  =  u1a/x1  is  the  Peclet  number  for  substrate  1.  At  high  Peclet  numbers  the 
aspect  ratio  for  conduction  is  large  and  conduction  in  the  x  direction  becomes  negligible  compared  with 
that  in  the  y  direction.  Heat  is  thus  transferred  into  the  body  by  conduction  and  carried  away  by  its  motion. 
The  same  argument  can  be  applied  to  body  2  if  it  is  moving  rapidly.  The  coating,  designated  as  region  3, 
will  have  a  high  aspect  ratio  since  it  is  assumed  throughout  the  analysis  that  a/h  >>  1.  The  energy 
equation  for  region  i  may  then  be  written  as: 


Piciui 


3T 

dX 
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(8) 


The  heat  generation  rates  per  unit  volume  R,  will  be  0  for  the  elastic  substrates  and  will  be  the  product  of 
the  shear  stress  and  the  inelastic  portion  of  the  shear  rate  for  the  coating.  Thus  from  Equations  (1)  and  (2) 
the  heat  generation  rate  in  the  film  is 


R3  -  { 


,  t'<x 

**  Tc 


T  *  T 


m 


(9) 


Continuity  of  flux  implies  continuity  of  K,3T/dy  at  y  =  0,  -h  for  0  <  x  <  2a  .  The  boundary  conditions  on 
the  temperature,  T,  are 
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lim  T  -  T0  and  T  -  T0  @  x  -  0  (10) 

|y|"“ 

Although  Equations  (8)  -  (10)  can  be  solved  in  their  present  form,  further  simplification  can  be  obtained  when 
the  thermal  resistance  of  the  coating  is  small  compared  with  that  of  either  of  the  substrates.  An  upper 
bound  to  the  temperature  rise  from  the  surface  to  the  center  of  the  film,  AT3,  may  be  obtained  by  neglecting 
the  effects  of  the  convective  motion  of  the  film.  The  temperature  rise  in  the  film  thus  corresponds  to  that 
due  to  conduction  between  two  parallel  plates  for  a  material  generating  heat  at  a  rate  R3  which  is 
AT3  =  R3h2/(8K3)  .  At  any  point  x  along  the  film,  the  heat  flux  going  into  region  1  is  q1  =  tyR3h,  where 
i|r  is  the  fraction  of  heat  generated  going  into  region  1.  One  may  obtain  an  estimate  of  the  surface 
temperature  rise  by  assuming  R3  and  hence  q,  to  be  independent  of  x.  The  thermal  problem  reduces  to 
that  for  a  semi-infinite  body  moving  past  a  uniform  heat  flux  starting  at  x  =  0  whose  solution  (or  that  for  the 
analogous  time  transient  problem  of  a  uniform  flux  applied  to  the  surface  of  a  stationary  semi-infinite  body 
at  time  t  =  o)  may  be  found  in  many  heat  transfer  texts  (ie.  Carsiaw  and  Jaeger3).  If  one  divides  AT3  by 
the  surface  temperature  rise  at  the  center  of  contact  (x  =  1)  obtained  from  the  solution  to  the  above 
problem,  AT^,  the  following  result  is  obtained: 

AT3  h  Kj  yjn  P 1 

at;  aK3  8n 

Although  the  Peclet  number,  P.,,  is  large,  h/a  will  be  very  small  for  thin  films  and  the  above  temperature  ratio 
will  be  small  as  long  as  the  thermal  conductivities  of  the  film  and  the  substrate  are  of  the  same  order  of 
magnitude.  This  is  frequently  true  for  solid  films,  especially  metallic  films  but  is  generally  not  true  for  liquid 
lubricants  which  have  much  lower  thermal  conductivities.  In  the  latter  case  temperature  variations  across 
the  film  can  be  significantly  larger  than  surface  temperature  rises. 

The  above  arguments  provide  the  criteria  for  neglecting  temperature  rises  in  the  film.  Doing  so,  reduces 
the  thermal  problem  from  a  three  region  to  a  one  region  problem  and  more  importantly  provides  justification 
for  property  variations  across  the  film.  In  the  remainder  of  this  analysis  temperature  rises  across  the  film 
will  be  neglected  but  the  ability  to  handle  the  three  region  thermal  problem  at  a  later  time  will  be  retained. 

With  the  above  approximations  Equation  (8)  needs  to  be  solved  only  in  region  1,  R,  =0,  and  the  boundary 
condition  at  y  -  given  in  Equation  (10)  is  replaced  by  the  condition  at  y  =  0  that  K.,aT/dy  =  h\|iR3. 
By  equating  temperatures  at  y  =  0  and  y  =  h  and  noting  that  q,  +  q2  =  hR3  ,  one  may  determine  the  flux 
partition  fraction 
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(11) 


One  may  now  introduce  the  reference  temperature  Tr  =  ursai|rTc0/K.,  and  the  additional  dimensionless 
variables  and  parameters 


f  -  —  ,  q,  -  — ^ —  ,  a  -  aph  ,  p  -  pph  ,  y  -  yTr  ,  X  -  XTr  .  (12) 

Tr  W*c0 


The  dimensionless  energy  equation  for  region  1  becomes 


at  _  i  a2f 

Pi  8r)2 


(13) 


with  the  boundary  conditions  on  t(5,ri) 


limt(5,ri)  «  t(0.ti)  -  0  (14) 


and  the  surface  flux  loading 


—  ?'<!>  (4-) 
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t|  “  0  ,  0  <  i  <  2  .  (15) 


C.  Solution  to  Combined  Problem 


Equations  (13)  to  (15)  provide  a  statement  of  the  thermal  problem  and  Equations  (4)  and  (5)  along  with  the 
prescribed  function  <t>,  represent  the  mechanical  problem.  They  are  coupled  through  the  property  variation 
relationships  given  in  Equation  (6)  which  may  be  expressed  in  dimensionless  form  as: 
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Q  _  &&P  -  yT  *  _  1  +  PP  ‘  _  v- 

^  ®  '  Tc  ~ =rr  ,  -  xxc 

1  +  AT 


where  p  =  p/ph  . 

A  direct  iteration  scheme  is  used  to  couple  the  thermal  problem  with  the  mechanical  problem.  Initially  the 
mechanical  problem  is  solved  for  isothermal  conditions  (t  =  0).  This  provides  an  initial  shear  stress 
distribution  i  and  from  Equation  (15)  an  initial  flux  loading,  qf(£) .  One  may  now  solve  the  thermal  problem 
for  a  new  temperature  distribution  to  evaluate  temperature  dependent  properties  from  Equation  (16)  which 
are  then  used  in  the  next  iteration  of  the  mechanical  problem.  Convergence  is  achieved  when  the  traction, 
defined  as  integral  of  the  shear  stress  across  the  contact  zone,  deviates  by  less  than  a  prescribed  tolerance 
between  iterations. 

At  each  iteration  Equation  (4)  may  be  treated  as  a  first  order  non-linear  differential  equation  below  the 
limiting  shear  stress  and  an  algebraic  equation  above  it.  At  low  pressures  the  Deborah  number,  D,  which 
multiplies  the  derivative,  can  become  quite  small  and  compared  with  the  second  term  on  the  left  side  of 
Equation  (4).  When  this  occurs  Equation  (4)  assumes  a  very  "stiff"  behavior  and  explicit  numerical 
integration  methods  such  as  Runge-Kutta  tend  to  become  unstable.  To  avoid  this  problem,  a  semi-implicit 
algorithm  was  implemented  that  has  been  widely  used  for  stiff  systems  (see  Press,  et  al.4)  and  is  described 
below. 

Put  the  top  part  of  Equation  (4)  in  the  form  di'/dt  =  F(Z,i') ,  discretize  the  (■  interval  of  interest  into 
a  grid,  denote  the  values  of  i  and  i'  at  the  jth  point  as  and  i'j  ,  and  use  the  average  value  of  F  to 
increment  i'  from  it  value  at  ^  to  that  at  5j+1  as  follows: 


V  *  F(5|,t'|)1  • 


,f  F^)+1>*j  +1j  in  the  above  equation  is  approximated  by 


F(5h,x'h)  -  F(5m,t,j)  +  (*',+1 


and  the  ensuing  expression  is  solved  for  t'j +1  ,  the  result  is 
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FKm.t7,)  +  Ffr.t'j) 


2 

3F 

dt7 

The  above  procedure  is  generally  stable,  provides  quadratic  accuracy  and  has  proven  to  be  very  effective 
for  the  problem  at  hand. 

Although  the  surface  temperature  T(£,0),  that  satisfies  Equations  (13)  and  (14),  for  any  surface  flux 
distribution  q.,(5).  can  he  determined  with  the  convolution  integral3 

t(5,0)  -  — !—  /-5^d{'  , 

o  vi7? 

it  was  found  more  expedient  to  adapt  an  existing  finite  difference  routine  written  by  the  author5.  The 
computer  time  is  approximately  the  same  and  the  finite  difference  approach  provides  sub-surface 
temperatures  without  the  need  for  further  computation.  In  addition  the  finite  difference  solution  is  readily 
amenable  for  use  in  treating  the  three  region  problem  discussed  previously  as  well  as  handling  more  general 
thermal  boundary  conditions  such  as  convection  or  radiation  from  the  surfaces. 

D.  Parametric  Studies 

A  number  of  parametric  studies  have  been  performed  to  show  the  influence  of  viscosity  and  shear  modulus 
on  the  nature  of  the  shear  stress  profiles  and  the  traction  vs.  slip  rate  behavior.  A  base  case  has  been 
selected  which  will  be  referred  to  as  the  low  viscosity,  low  compliance  case.  Medium  and  high  viscosity 
cases  will  correspond  to  a  factor  of  10  and  100  increase  in  viscosity  respectively.  Medium  and  high 
compliances  will  correspond  to  10  and  100  fold  decreases  in  shear  modulus.  Non  dimensional  parameters 
for  the  base  case  are  D  =  3.2x1  O'6  ,  S  =  6.4x1  O’3  s,  a  =  10  and  x  -  1  ■  The  above  parameters  would 
correspond  to  the  physical  parameters  ur  =  1000  in/sec,  h  =  6x10'5  in,  a  =  6x10"3  in,  tc0  =  6000  psi, 
G  =  1.2x105psi,  |i0  =  2.3x10"®  reyns,  a  =  10"4  in2/lb  and  ph  =  10s  psi. 

The  variation  of  the  average  shear  stress,  i,  with  the  slide  to  roll  ratio,  s,  is  shown  in  Figures  2  -  4  at  various 
compliances  for  low,  medium  and  high  viscosities  respectively.  (The  average  shear  stress  is  defined  as  the 
total  load  divided  by  the  total  contact  area  and  i  may  be  converted  to  a  friction  coefficient  by  dividing  it 


Figure  2.  Traction  vs.  slip  rate  at  low  viscosity, 
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model,  isothermal. 
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by  *ph/(4Tc0),  which  has  a  value  of  13.1  for  all  of  the  numerical  examples  presented  in  this  section.)  The 
curves  denoted  by  (1),  (2)  and  (3)  represent  low,  medium  and  high  compliance  respectively.  These  curves 
were  obtained  with  the  use  of  the  arctan  model  without  the  inclusion  of  thermal  effects.  The  arctan  model 
provides  an  asymptotic  approach  to  the  limiting  shear  stress  when  x  =  1  •  A  value,  of  x  =  -999  was 
selected  for  use  with  the  arctan  model  to  provide  a  smooth  transition  but  a  definite  cutoff  point  for  the 
limiting  shear  stress. 

Comparison  of  curves  (1)  and  (2)  in  Figure  2  reveals  that  an  order  of  magnitude  drop  in  shear  modulus 
produces  very  little  change  in  predicted  traction  at  low  viscosity  and  the  behavior  is  predominantly  visco¬ 
plastic  in  nature.  The  behavior  is  somewhat  different  at  high  viscosity.  Curves  (1)  and  (2)  in  Figure  3  differ 
significantly  at  low  slip  rates  but  approach  each  other  at  higher  slip  rates  where  plastic  behavior 
predominates.  Curve  (3)  in  Figure  4  lies  fairly  close  to  the  corresponding  curve  in  Figure  3  thus  indicating 
predominantly  elastic-plastic  behavior  where  an  order  of  magnitude  change  in  viscosity  has  very  little  effect 
on  traction.  At  high  viscosity,  low  shear  modulus  and  low  slip  rate,  the  behavior  becomes  primarily  elastic 
with  the  dimensionless  average  shear  stress,  f,  approaching  S/D  corresponding  to  integration  of 
Equation  (4)  with  <J>  set  equal  to  0. 

In  all  of  the  curves  described  above  the  average  shear  stress  increased  with  increasing  slide  to  roll  ratio. 
Inclusion  of  temperature  dependent  properties  in  some  circumstances  can  actually  show  a  decrease  in 
traction  due  to  a  drop  in  the  limiting  shear  stress  and  viscosity  as  the  coating  heats  up.  In  order  to  include 
thermal  effects,  the  additional  parameters  r|r.  Y.  X  and  P1  must  be  prescribed.  A  value  of  i|r  =  .5  was 
selected  for  the  thermal  partition  fraction  as  would  be  expected  at  low  values  of  s,  if  the  same  material  were 
used  for  both  substrates.  A  value  of  K,  =  6  lb-sec/0  F  corresponding  approximately  to  the  thermal 
conductivity  of  steel  was  selected  which  results  in  a  reference  temperature  of  Tr  =  3x1 03  s  °F  .  Values 
ofy  =  .00667  and  X  =  .0167  were  chosen  which  correspond  to  factor  of  two  drops  in  viscosity  and  limiting 
shear  stress  at  temperature  rises  of  100°F  and  60°  F  respectively  and  result  in  dimensionless  parameter 
values  of  y  =  20s  and  X  =  50s  .  Finally  a  value  of  the  Peclet  number  was  chosen  as  P,  =  100  for  the 
base  line  case  used  in  most  of  the  parametric  studies  with  a  value  of  P1  =  500  for  the  high  Peclet  number 
examples.  The  value  for  steel  would  lie  approximately  mid-way  between  these  two  cases.  The  effect  of  the 
Peclet  number  on  the  shear  stress  profiles  is  shown  in  Figure  5. 

The  variations  of  the  average  shear  stress  with  slide  to  roll  ratio  analogous  to  Figures  2  -  4  but  including 
thermal  effects  are  shown  in  Figures  6-8.  At  low  slip  rates,  thermal  effects  are  small  as  can  be  seen  by 
comparing  the  curves  in  Figures  6  -  8  with  the  corresponding  curves  in  Figures  2-4.  As  the  slip  rate 
increases  thermal  effects  cause  the  curves  to  level  off  more  rapidly  and  in  some  cases  tractions  are  seen 
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Figure  7.  Traction  vs.  slip  rate  at  medium  viscosity,  arctan  model  with  temperature  dependent 
properties. 


to  decrease  with  increasing  slip  rates.  At  higher  viscosities,  thermal  effects  become  visible  at  lower  slip 
rates. 

In  order  to  investigate  the  sensitivity  of  the  traction  behavior  to  the  rheological  model,  a  series  of  parametric 
studies  were  carried  out  with  the  use  of  the  Maxwell  model.  The  same  parameters  were  used  for  both 
models.  Comparisons  of  the  results  obtained  for  the  low  and  high  viscosity  cases  with  inclusion  of  thermal 
property  variations  are  shown  in  Figure  9  and  10.  The  curves  labelled  (1-a)  and  (3-a)  refer  to  the  arctan 
model  for  the  low  and  high  compliance  cases  respectively.  The  corresponding  curves  for  the  Maxwell  model 
are  labelled  (1-m)  and  (3-m).  The  Maxwell  model,  for  the  low  viscosity  case  is  seen  to  produce  slightly 
higher  tractions  since  the  creep  resistance  does  not  level  off  with  increasing  shear  rate  as  does  the  arctan 
model  (see  Figure  1 1 ).  For  the  high  viscosity  case  shown  in  Figure  10,  however,  the  behavior  is  very  nearly 
elastic-plastic  and  the  traction  curves  for  the  two  models  become  nearly  indistinguishable. 

The  arctan  model  avoids  the  abrupt  change  in  slope  of  the  stress  strain  rate  curve  at  the  yield  point 
characteristic  of  the  Maxwell  model.  This  abrupt  change  is  not  seen  in  the  traction  curves  shown  in  Figures 
9  and  10  because  it  is  smoothed  over  by  the  integration  process  used  in  evaluating  the  average  shear 
stress.  It  can  be  seen  clearly,  however,  in  the  shear  stress  profiles  shown  in  Figure  12.  The  upper  two 
curves  in  the  figure  do  not  include  the  temperature  variation  in  properties.  In  both  cases  the  abrupt 
transition  of  the  Maxwell  model  to  the  limiting  shear  stress  is  evident. 

The  corresponding  dimensionless  temperature  profiles  for  the  four  curves  in  Figure  12  are  shown  in  Figure 
13.  The  upper  curves  were  obtained  by  back  calculating  the  temperature  without  including  temperature 
variation  in  properties  which  corresponds  to  the  first  iteration  performed  in  the  process  of  obtaining  the 
thermal  solution.  At  the  slip  rate  of  s  =  .1  used  in  obtaining  Figures  12  and  13,  a  value  of  f  of  0.1 
corresponds  to  a  temperature  rise  of  30°  F. 
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III.  TRACTIONS  AT  HIGH  SLIP  RATES 

The  preceding  analysis  is  has  been  applied  to  situations  where  the  rolling  velocity  is  high  and  the  slip  rates 
are  relatively  low.  Under  these  conditions,  elasticity  effects  in  the  x  direction  can  be  important  and  a 
substantial  portion  of  the  shearing  can  occur  below  the  limiting  shear  stress.  Under  conditions  of  pure 
sliding  or  high  slide  to  roll  ratios,  yielding  will  occur  over  the  major  portion  of  the  contact  zone,  elastic 
strains  in  the  film  in  the  x  direction  will  be  relatively  small  and  visco-plastic  deformation  will  predominate. 
This  can  be  seen  at  slide  to  roll  ratios  as  low  as  5%  as  evidenced  by  the  close  proximity  of  curves  (1)  and 
(2)  in  Figures  2  -  4  under  a  factor  of  10  change  in  shear  modulus.  It  can  further  be  seen  from  Figure  4  that 
at  high  viscosity  the  dimensionless  average  shear  stress  at  a  10%  slide  to  roll  ratio  falls  within  5%  of  the 
value  of  t  =  1  which  would  occur  if  pure  plastic  behavior  (t  =  Tm)  were  to  prevail  over  the  entire  contact 
zone. 

A.  Equations  of  Plasticity  for  a  Thin  Film 

The  plastic  behavior  of  the  lubricant  may  be  analyzed  in  greater  depth  than  in  the  previous  sections  by 
treating  the  film  as  an  ideal  plastic-rigid  material  and  later  generalizing  the  analysis  to  include  visco-plastic 
behavior  as  used  previously.  The  Levy-Von  Mises  Equations  for  plane  strain  deformation  of  a  plastic-rigid 
material  will  be  presented  and  then  reduced  to  a  thin  film  form  based  an  a  small  aspect  ratio  (a/h0  <  <  1) 
analogous  to  the  reduction  of  the  hydrodynamic  equations  to  a  lubrication  environment. 

Development  of  the  Levy-Von  Mises  Equations  can  be  found  in  many  texts6,7  and  will  not  be  derived  here. 
They  consist  of  equilibrium  relationships,  stress  -  strain  rate  and  continuity  relationships  and  a  yield  criterion. 
They  will  be  presented  here  in  terms  of  stresses  and  instantaneous  velocities  under  conditions  of  plane 
strain. 

The  equilibrium  equations  are  given  below: 

+  ii  -  0  ,  (17) 

ax  ay 

ii  +  d-^L  -o  .  (18) 

ax  ay 

Equations  (17)  and  (18)  represent  equilibrium  in  the  x  and  y  directions  respectively. 
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An  effective  stress  or  flow  stress,  a,  is  used  in  both  the  stress  -  strain  rate  equations  and  the  yield  criterion. 
A  corresponding  effective  strain  rate,  e,  may  also  be  introduced  and  expressed  in  terms  of  the  instantaneous 
velocities.  These  are  both  defined  below: 


The  stress  -  strain  rate  equations  for  plane  strain  may  be  written  in  terms  of  the  velocities  as  follows: 


3u 

dx 


^(ox-oy) 


(21) 


dv 

ay 


6u  +  dv  _  2e 
ay  dx  a 


The  continuity  equation  may  be  obtained  by  adding  Equations  (21)  and  (22) 


(22) 

(23) 


du  |  dV 
dx  ay 

Finally  the  Von  Mises  yield  criterion  may  be  expressed  as 

a  -  v/3k  ,  (25) 

where  k  is  the  shear  yield  stress  for  the  coating. 

The  above  equations  will  now  be  reduced  to  a  thin  film  environment  by  estimating  orders  of  magnitudes  for 
their  various  terms  and  neglecting  terms  of  order  h0/a  compared  with  unity.  Since  the  quantities  a  and  h0 
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represent  reasonable  scale  factors  for  x  and  y  in  the  coating,  the  dimensionless  coordinates  5  and  n 
defined  previously  will  be  used  here.  A  characteristic  sliding  velocity  u0  will  be  introduced  which  could  be 
the  surface  velocity  for  pure  sliding  or  the  slip  velocity  for  a  rolling-sliding  contact  at  high  slide  to  roll  ratios. 
This  velocity  will  be  used  as  the  scale  factor  for  u.  An  appropriate  scale  factor  reflecting  the  order  of 
magnitude  of  v  would  be  u0h0/a  as  can  be  seen  from  inspection  of  the  continuity  equation.  The  reduced 
velocities  u  =  u/u0  and  v  =  av/(h0u0)  will  now  be  introduced  and  will  be  assumed,  along  with  the 
coordinates  5  and  t]  to  be  of  the  order  of  magnitude  of  1.  The  effective  strain  rate  §,  given  in  Equation  (20) 
may  be  written  in  terms  of  these  variables  as 


h0 


e 


u0 


/3  00 

2  0ri 


(26) 


The  right  hand  side  of  Equation  (26)  was  obtained  by  neglecting  terms  of  the  order  of  (h0/a)2  compared 
with  1. 


One  may  now  substitute  the  approximate  value  for  e  obtained  from  the  right  hand  side  of  Equation  (26) 
along  with  the  value  of  a  given  by  Equation  (25)  for  the  corresponding  quantities  in  Equation  (23)  and 
express  the  result  in  terms  of  the  reduced  variables  given  above  as  follows: 


00  +  0V  ^  00  T 
0t|  v  a  J  05  0r|  k 

Again  assuming  that  du/dn  is  of  order  1,  one  may  neglect  the  second  term  on  the  left  side  of  the  above 
equation  and  solve  the  resulting  expression  for  t  to  obtain  the  approximate  relationship 

x  -  k  ,  (27) 

where  the  departure  of  x/k  from  1  is  of  the  order  of  (h0/a)2.  Equation  (27)  is  extremely  important  in  that 
it  requires  the  coordinate  shear  stress  to  be  equal  to  the  shear  yield  stress  everywhere  in  the  coating  within 
the  contact  zone  where  plastic  deformation  occurs.  Some  of  the  ramifications  of  this  property  will  be 
discussed  later. 


A  reference  pressure  p0  may  be  used  as  a  scale  for  the  normal  stresses  which  can  be  the  maximum  Hertz 
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pressure,  ph  or  any  characteristic  pressure  associated  with  the  normal  loading  in  a  non-hertzian  contact. 
One  may  now  introduce  the  reduced  stresses  and  o„  in  Equation  (21)  and  perform  the  same 
substitutions  for  e  and  o  as  used  in  deriving  Equation  (27)  above  to  obtain  the  approximate  relationship 


80 


i*  8U 


(5-5v)^  ~4Sk^ 


an 


as 


where  ft  =  (h0/a)(k0/p0),  ICq  is  the  shear  yield  stress  of  the  coating  at  ambient  pressure  and  temperature 
to  be  used  as  a  characteristic  value  and  k  =  k/kQ. 


Since  the  reduced  stresses  and  partial  derivatives  in  the  above  equation  are  assumed  to  be  of  order  1  and 
the  coefficient  on  the  right  hand  side  is  small,  the  two  normal  stresses  are  approximately  equal  (  o  »  o  ). 
Furthermore,  since  for  the  conditions  of  plane  strain  assumed  here  oz  =  (ox  +  oy)/2,  a  nearly  hydrostatic 
state  of  stress  will  exist  where 


ox  -  oy  ~  oz  -  -p  .  (28) 

The  relative  error  in  Equation  (28)  will  be  of  the  order  of  fi*  =  max[(h0/a)2,6J.  All  of  the  arguments 
presented  thus  far  allow  k  to  be  a  variable  dependent  on  both  temperature  and  pressure. 

One  may  substitute  Equations  (27)  and  (28)  for  the  shear  stress  and  normal  stress  in  Equation  (18)  to  obtain 
the  relationship  for  the  pressure  variation  across  the  film 

-  5—  -  0  ,  (29) 

an  H 

where  i  =  z/k0  -  k.  The  above  equation  shows  the  usual  result  that  the  pressure  does  not  vary 
significantly  across  a  thin  film  along  with  an  indication  of  the  relative  error  (of  the  order  of  6*).  When  the 
same  substitutions  are  introduced  into  Equation  (17)  the  following  expression  is  obtained: 


dp  "  a  k  dx 
a$  h0  p0  dr) 


(30) 


Although  the  coefficient  of  di/dr]  will  be  large,  di/dr\  itself  will  be  very  small.  If  the  shear  yield  stress 
were  constant  then  k  =  1  and  from  Equation  (27),  i  =  1  +  terms  of  order  (h0/a)2  and  the  right  hand  side 
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of  Equation  (30)  would  be  of  the  order  of  6.  If  k  were  to  depend  on  pressure,  even  to  the  extent  that 
dk/3p  were  of  order  1,  the  pressure  variation  across  the  film  would  still  be  of  the  order  of  i*  and  the  right 
hand  side  of  Equation  (30)  would  at  the  most  be  of  the  order  of  max[(k0/p0)2,&]  which  will  still  be  small. 
The  above  arguments  indicate  that  under  ideal  plastic  conditions  the  pressure  in  the  contact  zone  for  a  thin 
film  under  high  shear  is  constant  or 


&  -  -  0  .  (31) 

as  an 

An  exception  to  the  above  relationship  would  occur  if  k  were  strongly  dependent  on  temperature  and  severe 
temperature  gradients  were  to  exist  across  the  film.  It  was  shown  in  the  first  section  of  this  report  that  such 
behavior  is  not  likely  to  occur  due  to  shearing  of  a  thin  metallic  film. 

B.  Deformation  of  an  Elastic  Substrate  Under  a  Uniform  Pressure 


The  requirement  that  the  pressure  must  be  constant  in  the  contact  zone  can  significantly  alter  the  shape  of 
the  contact  when  elastic  substrates  are  involved.  For  example,  based  on  plane  strain  Hertz  contact  theory, 
a  cylinder  of  radius  r  will  flattened  in  the  interval  0  <  x  <  2a  by  an  elliptical  pressure  distribution  having  an 
amplitude  ph  =  aE/[2r(1-v2)].  If  d  is  taken  to  be  the  displacement  of  the  surface  of  the  cylinder  relative 
to  a  line  drawn  tangent  to  the  deformed  cylinder  at  the  center  of  contact  (x=a),  then  d  will  be  zero  inside 
the  contact  zone  and  rise  sharply  outside  the  contact  zone.  A  dimensionless  displacement  defined  as 
d  =  2rd/a  is  plotted  against  £  for  both  the  hertzian  pressure  distribution  described  above  and  for  uniform 
pressures  over  the  interval  0  <  x  <  2a  in  Figure  14.  The  curve  labeled  (a)  represents  the  hertzian 
displacement.  The  curve  labeled  (b)  corresponds  to  the  displacement  of  the  same  surface  under  a  uniform 
pressure  that  integrates  to  the  same  load  as  the  hertzian  curve  (a).  The  curve  labeled  (c)  is  the 
displacement  at  the  uniform  pressure  required  to  flatten  the  cylinder  at  the  center  of  contact  which  may  be 
shown  to  be  twice  that  used  in  obtaining  curve  (b). 

Figure  1  depicts  a  situation  where  two  rollers  are  loaded  together  in  a  hertzian  contact  and  separated  by 
a  constant  coating  thickness  in  the  contact  zone.  This  is  generally  assumed  to  occur  when  the  coating 
thickness  is  smaller  than  the  overall  elastic  deformation  of  the  substrates.  The  shape  of  the  film  thickness 
profile  would  thus  be  the  same  as  the  sum  of  the  elastic  surface  displacement  profiles.  If  the  sliding  rate 
is  sufficiently  high  to  cause  a  significant  degree  of  yielding,  the  coating  will  not  support  a  hertzian  pressure 
profile  and  the  film  thickness  cannot  remain  constant.  The  variations  in  the  substrate  surface  displacements 
under  a  uniform  pressure  shown  in  Figure  14  can  be  large  compared  with  the  film  thickness  thereby 
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Figure  14.  Comparison  of  dimensionless  surface  displacement  profiles  under  elliptical  and  uniform 
pressure  distributions. 
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producing  major  film  thickness  variations. 

These  film  thickness  variations  have  not  been  treated  in  the  initial  traction  analysis  presented  in  Section  II 
or  in  the  elastohydrodynamic  traction  analyses  which  assume  limiting  shear  behavior  due  to  solidification 
to  occur  with  a  constant  film  thickness.  The  above  results  indicate  that  substantial  qualitative  behavioral 
discrepancies  may  result  from  separating  the  traction  theory  from  the  film  thickness  predictions  when 
solidification  or  limiting  shear  behavior  is  predicted  to  occur. 

C.  Example  of  a  Rigid  Slider  on  a  Substrate  with  a  Visco-Plastic  Coating 

Further  insight  into  the  effects  of  plastic  yielding  at  high  slip  rates  may  be  obtained  by  considering  a  visco¬ 
plastic  film  having  a  rheological  model  of  the  form 

x1  -  p—  ,  x  -  min(T/,k)  .  (32) 

dy 

A  schematic  of  the  film  under  consideration  is  shown  in  Figure  15.  if  a*  any  point  x,  the  shear  stress  in  the 
film  is  below  the  shear  yield  stress,  k,  then  the  equilibrium  equations  fc  the  two  dimensional  thin  film  will 
be  given  by 


iE  -  ii  ,  -  o  .  (33) 

dx  dy  dy 

If  the  film  is  assumed  to  be  steady  (w=0  at  y=h),  one  may  integrate  the  continuity  equation  given  by 
Equation  (28)  across  the  film  to  obtain  the  expression 

fudy  .  Q  -  ,  (34) 

{  2 

where  Q  denotes  the  volumetric  flow  rate  per  unit  transverse  length  and  h*  =  2Q/u0  . 

Since  from  Equation  (33)  the  shear  stress  must  vary  linearly  across  the  film,  it  will  have  its  maximum  and 
minimum  values  at  the  surfaces.  It  is  thus  possible  for  the  shear  stress  at  any  given  x  to  equal  the  shear 
yield  stress  on  one  of  the  surfaces  and  fall  below  the  shear  yield  stress  and  behave  in  a  viscous  fashion 
between  the  surfaces  thereby  supporting  a  pressure  gradient  along  the  film.  This  would  be  manifested  by 
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slip  at  the  yield  surface.  If  at  any  value  of  x,  the  shear  yield  is  reached  at  both  surfaces  (in  the  same 
direction)  then  the  shear  stress  will  be  constant  across  the  film  and  from  Equation  (33)  as  well  as  the 
arguments  presented  in  the  preceding  section,  the  pressure  gradient  along  the  film  must  vanish. 

For  the  conditions  of  pure  sliding  assumed  here,  the  boundary  conditions  on  the  velocity,  u,  under  usual 
conditions  of  no  slip  would  be  u=u0  at  y=0  and  u=0  at  y=h.  If  slip  were  to  occur  at  y=0  the  surface 
velocity  of  the  coating  may  be  represented  as  u0  -  Au0  similarly  the  velocity  may  be  represented  as  Au0 
for  slip  at  y=h. 

One  may  now  differentiate  Equation  (32)  with  respect  to  y,  substitute  Equation(33)  for  dr/dy,  integrate  the 
resulting  expression  subject  to  the  boundary  conditions  on  u  described  above  to  obtain  the  velocity 
distribution 


u  -  ~--^y(y-h)  +  (uo- Au0)i^  +  Auh-^  ,  (35) 

2p  dx  h  h 

which,  with  the  exception  of  the  slip  velocities,  is  part  of  the  procedure  generally  used  in  developing  the 
lubrication  equations. 


The  shear  stress  due  to  viscous  motion  may  be  obtained  by  differentiating  the  velocity  distribution.  The 
shear  stresses  at  y=0,  denoted  by  tq  and  at  y=h  by  rh  are  given  below: 


-To  - 


-Th  - 


h  dp 

2  ax 


h  dp 

2  ax 


(Au0  +  Auh) 

^  L. 


(Au0  +  Auh) 
h 


(36) 

(37) 


One  may  substitute  Equation  (35)  for  u  in  Equation  (34)  and  integrate  the  resulting  equation  across  the  film 
to  obtain  the  flow  relationship 


h3  ap 
i2|x  ax 


^0 

2 


(h-h*)  + 


(38) 


For  the  example  here  of  pure  sliding  with  a  monatonically  decreasing  film  thickness,  four  modes  of  shearing 
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will  exist  depending  on  the  values  of  the  surface  shear  stresses  -t0  and  -th  in  relation  to  the  shear  yield 
stress  k.  (It  should  be  noted  that  the  minus  signs  are  used  here  since  an  absolute  value  is  required  for 
comparing  the  shear  stress  to  the  shear  yield  stress  and  for  the  coordinate  system  in  this  example  with  the 
surface  at  y=0  in  motion,  the  shear  stresses  will  in  general  be  negative  when  Liieir  magnitudes  approach 
the  yield  value.)  These  will  be  described  below. 

Case  1.  No  slip,  -r0 ,  -rh  <  k.  The  motion  in  this  case  will  be  purely  viscous  with 
Au0  =  Auh  =  0  and  Equation  (38)  reduces  to 


I '  6"u° 


(h-h’) 

h3 


Case  2.  Slip  aty=0,  no  slip  aty=h,  -r0  =  k,  -rh  <  k.  One  may  set  -t0  =  k  and  Auh  =  0 
in  Equation  (36)  and  solve  for  Au0  to  obtain  relationship 


Au0  -  U0  + 


h2  dp 
2[i  dx 


kh 

n 


which  may  be  substituted  for  Au0  in  Equation  (38)  which  is  in  turn  solved  for  the  pressure 
gradient 


dp  _  3f*u0 
dx  2h3 

Case  3.  Slip  aty=h,  no  slip  aty=0,  -rh  =  k,  -r0  <  k.  One  may  set  -th  =  k  and  Au0  =  0 
in  Equation  (37)  which  maybe  solved  for  Auh  and  combined  with  Equation  (38)  as  in  Case 
2  to  obtain  the  following  relationships  for  the  slip  velocity  and  the  pressure  gradient: 


Auh 


U0  - 


h2  0p 
2p  dx 
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Case  4.  Slip  at  y=0  and  y=h,  -r0  =  k,  -rh  =  k.  Here  the  interior  of  the  film  will  be  plastic 
hence  the  pressure  gradient  must  vanish 


-  0  . 

3X 

The  relationships  for  the  pressure  gradients  and  surface  shear  stresses  for  the  various  cases  may  be  put 
in  dimensionless  form  with  the  use  of  the  following  additional  dimensionless  variables  and  parameters 


The  results  are  summarized  below: 


iP  . 
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A+h8p  Case1 

h  2  35 
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1  +  h-^  ,  Case  3 

as 

1  ,  Case  4 


i.fiap  Casel 

h  2  35 
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1  ,  Case  3 

1  ,  Case  4 


(39) 


(40) 


(41) 
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The  above  equations  provide  an  analytical  procedure  for  determining  the  pressure  gradient  and  surface 
shear  stress  which  is  combined  with  a  numerical  procedure  for  determining  the  value  of  h*  required  to 
satisfy  the  boundary  conditions  on  the  pressure  that  p  =  0  at  \  =  0,1  which  may  be  expressed 
mathematically  by  integrating  Equation  (39) 


f&di  -  0  .  (42) 

l  * 

The  solution  algorithm  first  involves  dividing  the  interval  $1  into  a  discreet  number  of  points  and 
selecting  a  trial  value  of  h*.  Starting  at  the  first  point  one  may  now  compute  the  pressure  gradient  and 
corresponding  surface  shear  stresses  using  the  Case  1  form  of  Equations  (39)  -  (41)  denoting  the  maximum 
of  i0  and  Th  as  xmax  for  Case  1.  If  rmax  s  1,  one  may  save  the  value  of  the  pressure  gradient  and 
proceed  to  the  next  point.  If  for  Case  1,  Tmax  >  1  and  3p/3$  >  0,  the  solution  will  either  correspond  to 
Case  2  or  Case  4.  To  determine  which,  compute  the  pressure  gradient  and  shear  stresses  for  Case  2.  If 
ih<;1  proceed  to  the  next  point.  If  not,  set  the  pressure  gradient  and  shear  stresses  to  their  Case  4  values 
and  proceed  to  the  next  point.  If  for  Case  1  tmax  >  1  and  dp/31;  <  0  the  solution  will  either  correspond 
to  Case  3  or  Case  4.  To  determine  which,  compute  the  pressure  gradient  and  shear  stresses  for  Case  3. 
If  t0  s  1  proceed  to  the  next  point.  If  not,  set  the  pressure  gradient  and  shear  stresses  to  their  Case  4 
values  and  proceed  to  the  next  point. 

The  above  procedure  may  be  used  to  determine  the  pressure  gradient  at  every  point  which  may  then  be 
used  for  numerical  integration  of  the  left  side  Equation  (42).  The  correct  value  of  h*  will  be  the  one  that 
satisfies  Equation  (42).  For  the  present  analysis,  numerical  integration  was  carried  out  using  the  trapezoidal 
rule  and  h*  was  determined  using  the  Brent  method3.  The  combined  analytical  and  numerical  procedure 
used  here  is  a  reduction  of  a  general  purpose  numerical  algorithm  developed  by  the  author  for  treating 
classes  of  one  dimensional  non-newtonian  rolling-sliding  contacts  in  powder  lubrication8. 

Once  the  pressure  gradients  are  determined  they  may  be  integrated  with  respect  to  5  to  compute  the 
pressure  distribution  which  may  in  turn  be  integrated  to  predict  the  dimensionless  load  W.  One  may  also 
integrate  t0  given  by  Equation  (40)  to  compute  the  tractive  force.  These  dimensionless  forces  are  defined 
as 
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A 

w  - 


ho 

a2k 


where  W  and  f  denote  normal  and  frictional  forces  per  unit  transverse  length. 


Computations  have  been  carried  out  for  the  case  of  a  convergent  slider  with  a  linear  film  thickness  profile 
given  by 


h  -  h,  -  (h,  - 1)5  ,  0*5  si  . 

The  variation  of  the  dimensionless  load  capacity,  W  and  traction,  f,  with  the  viscosity-speed  parameter  4. 
are  given  in  Figures  16  and  17  respectively  at  various  values  of  the  film  thickness  ratio  h,.  The  parameter 
4  is  an  index  of  the  relative  magnitude  of  the  viscous  shear  force  to  the  yield  stress  and  is  proportional  to 
both  the  viscosity  and  the  sliding  speed.  At  low  values  of  (1  the  flow  will  be  entirely  viscous  in  nature  with 
shear  stresses  below  the  yield  stress  everywhere  in  the  film.  At  the  points  denoted  by  the  small  dots  (•) 
in  Figure  16,  the  shear  yield  stress  is  reached  on  the  stationary  surface  at  the  trailing  edge  where  the  film 
thickness  is  smallest.  As  4  is  increased  further  the  general  level  of  shear  goes  up  in  the  film  and  the  zone 
at  which  trailing  edge  slip  occurs  increases.  At  the  points  denoted  by  the  circles  (°)  the  shear  yield  stress 
is  reached  and  slip  starts  to  occur  on  the  moving  surface  near  the  leading  edge  of  the  film.  As  4  continues 
to  increase,  both  slip  zones  expand  inward  until  they  eventually  overlap.  Points  of  incipient  overlap  are 
indicated  by  the  intersections  of  the  various  curves  with  the  broken  line  in  Figure  16. 

The  behavior  of  the  film  at  values  of  4  to  the  right  of  the  broken  line  can  be  observed  by  examining  the 
pressure  profiles  shown  in  Figure  18  for  the  particular  case  of  h1  =  3  .  The  top  curve  corresponds  to  the 
threshold  value  of  4  where  yield  across  the  film  occurs  at  one  point.  As  (1  increases  the  fully  plastic  zone, 
as  manifested  by  the  flat  portions  of  the  curves,  expand  and  the  pressures  drop.  As  the  sliding  speed  (or 
viscosity)  becomes  high  enough  the  plastic  zone  extends  to  the  entire  film  and  the  pressure  (hence  the  load 
capacity)  falls  to  0.  As  can  be  seen  from  Figure  16,  this  complete  loss  of  load  capacity  occurs  at  finite 
defined  values  of  4  which  depend  only  on  the  film  thickness  ratio.  This  value  of  4  which  will  be  denoted 
by  4*  may  be  computed  from  Equation  (39)  for  any  smoothly  convergent  film  as  shown  below. 

Just  before  the  entire  film  becomes  plastic  a  small  region  will  exist  in  the  inlet  where  Case  2  shear  prevails 
and  a  region  at  the  exit  where  Case  3  shear  prevails.  The  pressure  gradient  must  approach  0  as  each  of 
these  regions  merge  with  the  fully  plastic  region.  As  the  length  of  the  Case  2  inlet  region  becomes  very 
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small  the  film  thickness  h  will  approach  hv  By  replacing  |i  with  |im,  h  with  hv  and  setting  the  pressure 
gradient  to  0  for  Case  2  of  Equation  (39)  one  obtains  the  relationship 


(43) 


One  may  apply  a  similar  argument  to  the  exit  region  and  replace  ji  with  jim,  h  with  1,  and  set  the  pressure 
gradient  to  0  for  Case  3  of  Equation  (39),  then  use  Equation  (43)  to  eliminate  h*  and  solve  the  resulting 
expression  for 


A* 


1  +fi? 
2 


(44) 


Thus  for  any  given  inclination,  Ah  =  h1  -  h0  and  minimum  film  thickness  h0,  the  load  capacity  will  vanish 
when  the  product  of  the  viscosity  and  sliding  speed  become  high  enough  so  that  £  ;>  ji*  The  variation  of 
film  thickness  with  speed  or  viscosity  at  fixed  load  and  inclination  can  be  obtained  by  plotting  the  variation 
of  h0/Ah  =  1/(h.,-l)  with  jiu0/(kAh)  =  p./(h1-1)  at  fixed  values  of  WAh/(a2k)  =  W(h.,-1).  This  may  be 
obtained  from  cross  plots  of  the  type  of  data  used  in  generating  Figure  16  and  is  shown  for  the  particular 
case  of  W(ht-1)  =  ,04,  in  Figure  19. 

The  film  thickness  as  shown  by  solid  curve  in  Figure  19  is  seen  to  increase  with  increasing  speed  (or 
viscosity)  at  low  speeds  then  reach  a  maximum  and  fall  off  suddenly  as  plastic  yielding  starts  to  occur.  The 
minimum  film  thickness  then  continues  to  decrease  but  at  a  less  drastic  rate.  The  ratio  of  the  inclination  of 
the  trailing  edge  film  thickness  at  the  right  end  of  the  solid  curve  is  16  to  1  which  is  essentially  an  edge 
loading  condition.  In  a  separate  computation  stiffnesses  were  calculated  and  found  to  be  negative  for  the 
portion  of  the  solid  curve  having  a  positive  slope  and  lying  to  the  right  of  the  maximum,  indicating  an 
unstable  region.  In  a  transient  analysis  one  would  expect  a  rapid  drop  in  film  thickness  in  this  region  at  a 
rate  controlled  by  squeeze  film  effects.  The  dashed  line  in  the  figure  indicates  the  film  thickness  vs.  speed 
behavior  that  would  be  predicted  to  occur  in  the  absence  of  plastic  yielding  and  demonstrates  the  drastic 
effect  that  plastic  yielding  is  predicted  to  have  on  film  thickness  at  high  viscosity  or  speed. 
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IV.  SUMMARY  AND  CONCLUSIONS 

Analyses  are  developed  to  predict  the  tractive  forces  that  will  occur  in  a  rolling-sliding  contact  between  one 
body  of  a  single  material  and  a  second  substrate  having  a  thin  coating  of  solid  lubricant.  The  materials  are 
assumed  to  be  isotropic  and  conditions  of  plane  strain  are  assumed  to  prevail.  The  film  is  taken  to  be  thick 
compared  with  asperity  heights  but  very  thin  compared  with  contact  lengths.  The  coating  is  assumed  to 
be  soft  compared  with  the  substrates  and  the  contact  pressures  are  assumed  to  be  high  compared  with  the 
shear  yield  stress  of  the  coating. 

An  initial  traction  analysis  is  given  in  Section  II  which  includes  elastic,  plastic,  creep  and  thermal  effects. 
The  analysis  was  designed  to  treat  rolling  dominated  rolling-sliding  contacts  and  incorporates  an  approach 
developed  for  treating  solidification  in  elasto-hydrodynamic  contacts  but  oriented  toward  the  characteristics 
of  solid  films  with  the  inclusion  of  a  suitable  thermal  model.  Parametric  studies  have  been  performed 
showing  the  interaction  of  the  various  effects  for  two  rheological  models.  The  analysis  has  additional 
flexibility  that  has  not  yet  been  exploited  such  as  capability  of  handling  more  general  thermal  boundary 
conditions,  inclusion  of  non-hertzian  pressure  profiles  and  other  rheological  models.  It  is  also  highly 
tractable  and  has  the  capability  being  extended  to  handle  real  three  dimensional  contacts  with  variable  slip 
rate  etc.,  characteristic  of  real  machinery. 

An  analysis  was  formulated  for  a  coating  under  sliding  dominated  conditions  in  Section  III.  Elasticity  effects 
for  the  coating  were  neglected  and  the  Levy-Von  Mises  Equations  for  a  plastic-rigid  system  were  reduced 
based  on  thin  film  approximations.  The  analysis  shows  that  the  coating  will  both  be  in  a  hydrostatic  state 
of  stress  and  will  not  in  general,  be  able  to  support  a  pressure  gradient  once  yielding  occurs.  This  latter 
result  imposes  serious  restrictions  on  the  assumption  of  a  constant  film  thickness  used  in  Section  I,  as  well 
as  References  1  and  2  and  related  elastohydrodynamic  traction  analyses  that  use  limiting  shear  theory,  at 
high  slip  rates. 

A  simple  viscoplastic  model  was  applied  to  a  coating  on  a  rigid  substrate  moving  under  a  converging  rigid 
slider  which  permitted  yielding  within  the  film  or  at  the  coating-substrate  interfaces  together  with  viscous 
motion  of  the  coating  at  stress  levels  below  the  yield  point.  When  yielding  occurred  at  both  surfaces  at  any 
given  location  then  the  film  would  be  at  a  completely  plastic  state  at  that  location  hence  unable  to  support 
a  pressure  gradient.  The  results  show  that  for  any  given  lubricant,  inclination  and  minimum  film  thickness 
there  is  a  critical  speed  that  when  exceeded  results  in  a  complete  loss  of  load  capacity.  At  a  fixed  load  and 
inclination,  the  film  thickness  is  shown  to  drop  off  rapidly  as  significant  yielding  occurs  resulting  in  extremely 
small  film  thicknesses  compared  with  those  predicted  in  the  absence  of  yielding. 
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In  the  case  of  hertzian  contact,  if  the  viscosity  is  highly  pressure  dependent  as  assumed  in  Section  II  and 
the  shear  rates  are  sufficiently  low  so  that  yielding  occurs  over  a  relatively  small  portion  of  the  contact  zone, 
near  hertzian  pressures  (exception  for  a  small  amount  of  flattening  near  the  center  of  contact)  hence  a 
relatively  constant  film  thickness  in  the  contact  zone  can  be  maintained.  As  the  shear  rate  becomes 
sufficiency  high  to  cause  yielding  to  occur  over  a  major  portion  of  the  contact  zone  a  large  constant 
pressure  region  will  develop  and  sharp  variations  in  the  film  will  occur  in  the  case  of  compliant  substrates 
of  thin  films  which  make  it  inappropriate  to  separate  traction  and  film  thickness  theories. 
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V.  RECOMMENDATIONS 

The  conclusions  of  the  work  performed  thus  far  indicate  that  it  is  possible  to  predict  tractions  including 
elastic,  plastic,  creep  and  thermal  interactions  under  suitable  conditions  but  more  work  is  needed  to  extend 
the  results  further  into  the  plastic  regime.  In  particular  a  combined  film  thickness  an<_  traction  theory  is 
needed  to  extend  the  present  analysis  as  well  as  previous  elastohydrodynamic  analyses  to  account  for  the 
effects  of  lubricant  plasticity  on  the  film  thickness  profile. 

While  the  results  of  this  study  indicate  the  feasibility  of  a  continuum  approach  to  certain  modes  of  solid 
lubrication,  the  scope  needs  to  be  enlarged  to  treat  non-isotropic  lubricants  such  as  MoS2,  as  well  as  three 
dimensional  and  transient  contacts  characteristic  of  real  machinery. 

Predictions  of  small  film  thicknesses  do  not  necessarily  imply  failure  of  the  solid  lubricant.  When  the  local 
coating  thickness  is  of  the  same  order  of  magnitude  or  smaller  than  asperity  heights  on  the  substrate,  they 
will  interact  both  with  the  coating  and  with  each  other.  These  interactions  need  to  be  considered  in 
conjunction  with  the  lubricant  model  to  provide  a  general  theory  that  bridges  the  gap  between  thick  films 
and  classical  friction  theory. 

Experimental  verification  is  badly  needed  both  regarding  the  predicted  variations  in  film  thickness  at  high 
slip  rates  and  the  range  of  validity  assumptions  used  in  developing  the  initial  traction  model  in  Section  II. 
It  will  also  provide  guidance  for  treating  the  recommended  theoretical  extensions.  In  addition,  the  various 
lubricant  models  require  property  data  which  may  be  extracted  from  experiments  under  laboratory 
conditions  simulating  actual  contacts  when  primary  laboratory  data  are  not  readily  obtainable. 
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